"Chemical" composition of the Quark-Gluon Plasma in relativistic heavy-ion collisions 
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We study the evolution of the quark-gluon composition of the plasma created in ultra- Relativistic 
Heavy Ion Collisions (uRHIC's) employing a partonic transport theory that includes both elastic 
and inelastic collisions plus a mean fields dynamics associated to the widely used quasi-particle 
model. The latter, able to describe lattice QCD thermodynamics, implies a "chemical" equilibrium 
ratio between quarks and gluons strongly increasing as T — > T c , the phase transition temperature. 
Accordingly we see in realistic simulations of uRHIC's a rapid evolution from a gluon dominated 
initial state to a quark dominated plasma close to T c . The quark to gluon ratio can be modified by 
about a factor of ~ 20 in the bulk of the system and appears to be large also in the high pr region. 

We discuss how this aspect, often overflown, can be essential for a quantitative study of several 
key issues in the QGP physics: shear viscosity, jet quenching, quarkonia suppression. Furthemore 
a bulk plasma made by more than 80% of quarks plus antiquarks provides a theoretical basis for 
hadronization via quark coalescence. 

PACS numbers: 12.38. Mh, 25.75.Nq 



The search for the Quark-Gluon Plasma (QGP) started 
its golden age thanks to the experiments at the Relativis- 
tic Heavy Ion Collider (RHIC) that have supplied con- 
vincing physical observations that a new state of matter 
has been created [1, 2]. Such a matter has a very small 
viscosity [3-5], a high opacity to high— pr particles [6], 
a strong screening of the interaction able to significantly 
dissociate charmonia [7], and exihibits a modification of 
the hadronization respect to the vacuum toward a quark 
coalescence mechanism [8-10]. Furthermore some RHIC 
data hints to the creation of an "exotic" initial state of 
matter, the Color Glass Condensate, that could be the 
high-energy limiting state of the QCD interaction [6]. 
The new and upcoming experiments at the Large Hadron 
Collider (LHC) have confirmed the main gross properties 
observed at RHIC [11], but they will allow to explore a 
larger temperature range also with quite different heavy- 
quark abundancies and will provide more suitable condi- 
tions for creating CGC phase as initial state. 

The several probes mentioned above rely on the com- 
parison between experimental data and model predic- 
tions. A closer look into the several theorethical ap- 
proaches to the different QGP probes reveals that in 
some cases the QGP is described as a Gluon Plasma be- 
cause initially this should be the case. In fact most of 
the particles come from low x momentum fraction where 
the nucleon parton distribution functions are gluon dom- 
inated. This for example is the case of the most popular 
jet quenching models assuming a bulk gluon matter. In 
other cases as in the viscous hydrodynamics a chemical 
quark-to-gluon equilibration is implicit in the employ- 
ment of a lattice QCD (1QCD) equation of state P(e). 
For the study of quarkonia instead usually one considers 
a plasma of gluons or an equilibrated QGP according to a 
massless quark-gluon description acting for dissociation. 
On the contrary the observation of quark-number scaling 
in the elliptic flow and the baryon over meson enhance- 



ment are explained by quark coalescence models based 
on a quark dominance in the plasma [8, 10]. 

Certainly despite a lack of full integration of the dif- 
ferent descriptions of the QGP, all of these have been 
useful simplifications that allowed to successfully identify 
the creation of the QGP plasma and its gross properties. 
Nonetheless once we have identified the main qualitative 
features of the matter created in uRHIC's a quantita- 
tive knowledge of properties like the rj/s or the solution 
to open issues on the jet quenching mechanism, quarko- 
nia dissociation-regeneration, hadronization mechanism 
and existence of a CGC matter, requires to consider the 
poorly explored issue of the "chemical" composition of 
the QGP. 

The assumption of chemical equilibrium of the QGP, 
when considered, is usually discussed treating the QGP 
as a gas of massless quarks and gluons; therefore the ex- 
pected ratio is given simply by the ratio of the degrees of 
freedom N q+1 /N g = d q+ - q jd g = 2N c N f /{N 2 c - 1) = 9/4 
for a system with 3 flavors. On the other hand, as well 
known from 1QCD, the QGP appears to be significantly 
different from a mere massless gas, showing deviation 
of both the energy density and the pressure from the 
e/T 4 = 3P/T 4 = cost., and in particular exihibiting a 
large trace anomaly (9^) = e — 3P. It has been shown 
that such a behavior can be described in terms of a mas- 
sive quasiparticlc model in which both gluons and quarks 
acquire a thermal mass m(T) ~ g(T) T, as suggested 
also by the Hard Thermal Loop (HTL) approach [12-16] 
or dimensionally reduced screened perturbation theory 
(DRSPT) [17, 18]. 

In particular a quasi-particle model (QPM) with a 
T— dependent Bag constant has been successfully applied 
to quantitatively describe 1QCD results for equilibrium 
thermodynamics[19-21] including the recent ones per- 
formed with an unprecedented level of accuracy at the 
physical quark masses [22]. We notice that if the QGP 
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can be described in terms of finite mass excitations this 
has a strong impact on the quasiparticle chemical ratio 
N q+ q/N g . In fact at equilibrium one has: 

N q+q = d q+q m 2 q (T)K 2 (m q /T) 
N g d g m 2 (T)K 2 (m g /T)' ( > 

where K 2 is the Bessel function and m 9jff (T) are the 
T— dependent quark and gluon masses that can be de- 
termined by a fit [19] to recent 1QCD calculations [22]. 
In Fig.l, we show by solid line the equilibrium ratio 
when the fit to 1QCD e(T) is done assuming m 2 /m 2 = 
3/2 • {N 2 - l)/N c (2N c + N f ) = 4/9 according to a pQCD 
scheme [20]. We plot also by dash-dotted line the ex- 
pected ratio assuming the HTL m 2 /m 2 = 1/9 ratio. This 
of course leads to a larger quark abundancy. Even if in 
Ref. [19] some of the authors have shown that in the 
last case one can better describe the diagonal quark sus- 
ceptibilities, in the following we will consider the more 
commonly assumed pQCD case, which also prevents from 
overestimates of the magnitude of the dynamics studied. 
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FIG. 1: Quark to gluon ratio at equilibrium as a function 
of the temperature as predicted by QPM [19]. Solid line is 
assuming a m q /m g ratio according to pQCD and dot dashed 
line according to HTL; by dashed line it is indicated the mass- 
less case. 

In this Letter we discuss the issue of the quark-to-gluon 
ratio of the matter created in uRHIC's at both RHIC and 
LHC energies. We employ a Boltzmann-Vlasov transport 
theory to simulate the partonic stage of the HIC in a re- 
alistic way. In the last years several codes have been 
developed based on transport theory at the cascade level 
[23-25, 27], i.e. including only collisions between mass- 
less partons, with quite rare exceptions [28, 29]. These 
approaches have been more recently developed to study 
the dynamics of the partonic stage of the HIC at fixed 
shear viscosity [27, 29-31] with the advantages to explore 
possible effects of kinetic non-equilibrium, having also a 
wider range of validity in px and in rj/s. 

We present here for the first time results within a trans- 
port approach that includes the mean field dynamics as- 
sociated to the thermal self-energies generating the finite 



mass m(T) in the quasi-particle model discussed above 
[19, 20]. The approach is formally similar to the one de- 
veloped in Ref. [29] for the NJL mean field dynamics, but 
here the quasiparticle mean field allows to account for the 
proper equation of state, P(e), as evaluated from 1QCD. 
In such a picture the relativistic Boltzmann-Vlasov equa- 
tion can be written as follows: 

tp"0„ + m*{x)d^m*{x)d$]f{x,p) = C[f](x,p) (2) 

where C(x,p) is the Boltzmann-like collision integral, 
main ingredient of the several cascade codes: 

C=[f ftfvfo - hf2)\Mv2^12\ 2 5\ P l +P2-P 1 -P 2 ) 
2 1'2' 

(3) 

where J\ = J. d 3 pj /(2w) 3 2Ej, fj arc the particle distri- 
bution functions, while M. denotes the invariant tran- 
sition matrix for elastic as well as inelastic processes. The 
elastic processes have been implemented and discussed 
in several previous works [23-25, 27]. The inelastic pro- 
cesses between quarks and gluons {gg O qq) is instead 
the main focus of the present Letter. We have evaluated 
the matrix element in a pQCD LO order scheme. The 
tree diagrams contributing to the gg O qq correspond 
to the u, t, s— channels: M = M s + M t + M u . For 
the massless case the cross sections for such processes 
are the textbook pQCD cross section for jet production 
in high-energy proton-proton collisions. With massive 
quarks the calculations are the renowed Combridge cross 
sections [32] used to evaluate heavy quark production. 
In our case we have considered a finite mass for both 
gluons and quarks together with a dressed gluon prop- 
agator. The details of the calculations are lengthy and 
will be published elsewhere [34], and are quite similar to 
the one in [33] for finite current strange quark mass, we 
only mention that the cross section is dominated by the 
t— and u— channel while the s— channel and all the inter- 
ference terms are neglible. The squared matrix element 
of the t— channel is given by: 

\(t-m 2 q )M t \ 2 = ^ir 2 a 2 s [{m 2 q -t)(m 2 q -u) 

~2m 2 q (t + m 2 q )-4m 2 q m 2 g -m g }(4) 

and of course by crossing symmetry the u— channel can 
be obtained by u -H> t exchange. 

At the temperature explored by the system created 
at RHIC and LHC the energy scale of the processes, 
Q 2 = (2ttT) 2 , is of the order of - 4 - 6 GeV 2 there- 
fore the perturbative estimate of the strong coupling is 
a s ~ 0.7 — 0.8. Nonetheless the pQCD scheme with 
renormalized fermionic and bosonic lines should provide 
a useful guideline for the estimate of the pertinent cross 
section, keeping in mind that it can be expected that 
the real cross section could be even larger than the esti- 
mated ones. This is known to be the case for the elas- 
tic scattering processes. However this can only increase 
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the magnitude of the effect that we discuss here. Fur- 
thermore the evaluated cross sections are already large 
enough that their variations within a factor of two only 
marginally affects the final quark to gluon ratio, in fact 
the system gets anyway quite close to equilibrium. 

The thermodynamical self-consistency of the QPM re- 
quires a sclf-consitency between the Bag constant and 
the effective mass of the quasiparticles [19] which leads 
to a gap- like equation coupled to Eq.(2): 

dB _ r d 3 f mi (x) 

dm,' J (2tt)3 Ek{x) MX ' P) (5) 

with i = q,q, g. Eq.(5) allows to evaluate locally the mass 
in Eq.(2) also in non-equilibrium conditions guaranteeing 
the conservation of the energy-momentum tensor of the 
fluid. For the numerical solutions of Eqs.(2) and (5) we 
use a three dimensional lattice that discretizes the space 
and the standard test particle methods that samples the 
distribution function f(x,p) by means of an ensemble 
of points in the phase-space, for more details see Ref. 
[29]. We have carefully checked that the numerical im- 
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FIG. 2: Quark to gluon ratio as a function of time (normalized 
to r ) in Au + Au at y/s = 200AGeV (black lines) and for 
Pb + Pb at y/s = 5.5ATeV (light lines). Dashed lines are for 
the massless case and the solid for the massive case 

plementation in a stationary box is able to reproduce 
correctly both the kinetic equilibrium and in particular 
the chemical one, i.e. the right abundancy of quarks, an- 
tiquarks and gluons according to Eq.(l), which implies 
also a proper implementation of the detailed balance. 

We have simulated Au + Au collisions at y/s = 200 
AGeV and Pb + Pb at y/s = 5.5ATcV for 0-10% central- 
ity. The initial conditions in the r-space are given by the 
standard Glauber condition while in the p-space we use 
a Boltzmann-Juttner distribution function up to a trans- 
verse momentum px = 2 GeV while at larger momenta 
minijet distributions are implemented, as calculated by 
pQCD at NLO order [8]. At RHIC the initial maximum 
temperature at the centre of the fireball is Tq = 340 MeV 
and the initial time tq = 0.6 fm/c (corresponding also to 



the To • To ~ 1 criterium) as in succesfull hydrodynamical 
simulations. For Pb + Pb collisions at y/s — 5.5 ATeV we 
have T = 600 MeV and t ~ 1/T = 0.3 fm/c. In Fig.2 it 
is shown the time evolution of the ratio R qg = N q+q /N g 
for both the massless case (dashed line) and the massive 
quasi- particle case (solid line). We clearly see that in 
the massless case the system reaches very quickly, in less 
than 1 fm/c, the chemical equilibrium given by R qg ~ 2. 
Therefore in the massless case to assume a chemical equi- 
librium in modeling the plasma composition can be al- 
ready considered a reasonable approximation. 

When the quasi-particle massive case is considered, 
we can again see that one can reach quickly a value of 
R q g ~ 2.5, then there is a slower rise that leads to a 
continue increase up to R qg ~ 4.5 for Au + Au . We 
notice that for the massive case the equilibrium value is 
strongly T dependent, especially close to T c (see Fig.l), 
and the system more dilute in this stage is not capable 
to follow such a rapid change. Nonetheless we find that 
the fireball reaches a value relatively close to the equilib- 
rium at T ~ T c and eventually it is composed by about 
85% of quark plus anti-quarks. At LHC the trends are 
very similar but a longer part of the lifetime is spent in 
a T-region where the equilibrium R qg is nearly constant. 
This results into a moderately smaller final ratio. 

We show in Fig.2 the result starting from an initial 
gluon dominated plasma with a R qg =0.25 which is al- 
ready a quite strongly gluon dominated plasma, however 
changing this initial ratio by a factor of two affects the 
final ratio by less than a 10%, while the effect we are de- 
scribing modifies the value of R qg by more than an order 
of magnitude. 

In Fig. 3 the px dependence of the quark to gluon ratio 
is shown for the initial distribution (dotted line) and the 
freeze-out distributions: massless case (dashed line) and 
massive case (solid line). Black lines are for Au + Au 
and light ones fore Pb + Pb. We see the large differ- 
ence between the massless and the massive case and also 
that the net gluon to quark conversion extends up to 
quite large pr- The decrease of the ratio with pt can be 
expected considering that high— px particles can more 
easily elude the equilibration dynamics. However, in the 
massive case, we note a quite strong dependence below 
Pt = 2 GeV that has not to be interpreted as a fast 
detachment from the chemical equilibrium. In fact the 
dependence of the ratio can be evaluated analitically 
at equilibrium and it is given by 

dN/d 2 p T \ q+q = d q+ - q ?7 4 e 7[("4-^WT)/T] 
dN/d 2 p T \ g d g m 9 T e^ m T-PoVT)/T] <■ > 

where /3q is the radial flow velocity, m T and m T are 
the transverse masses. We plot in Fig. 3 by thin solid line 
such a function rescaled by a factor 0.85 accounting for 
the lack of full thermalization. The strong px dependence 
obtained in the transport simulation follows very closely 
the equilibrium behavior at least up to pt ~ 1.5 GeV. 
This is a well known effect predicted by hydrodynamics 
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and experimentally observed from SPS to LHC energy 
for hadronic spectra. 
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FIG. 3: Quark to gluon ratio vs pr at the freeze-out time. 
Labels as in Fig. 2. The thin solid line represents the full 
equilibrium raio, Eq. (6); see text for details. 



happens also at pr > 4 — 5 GeV in the region where 
the dominant hadronization mechanism should be inde- 
pendent fragmentation. In this region even if far from 
chemical equilibrium we still find a modification of R qg 
of about a factor of 7 respect to modelings discarding 
the q f-> g conversion mechanism. In Fig. 4 we show how 
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Our study shows that the most common microscopic 
description of the 1QCD thermodynamics implies a 
plasma that in uRHIC's evolves from a Glasma toward 
a quark dominated plasma. The R qg ratio can evolve by 
more than a factor 20 and at freeze-out is anyway more 
than a factor 2 larger than the one for an equilibrated 
massless QGP. Such an effect appears to be quite solid 
despite uncertainties in the inelastic cross sections and 
could be even larger and robust if one considers that the 
quark susceptibilities computed on 1QCD hints at even 
smaller m q /m g ratio [19] and hence quite larger equi- 
librium value of Rqg, see dashed-dotted line in Fig.l. 
Therefore the evolution of quark composition of the QGP 
should not be discarded for any quantitative studies of 
the several probes of the QGP and its properties. 

Our result provides also a support to the quark coales- 
cence hadronization mechanism, capable to explain two 
main observations at RHIC: the baryon/meson anomaly 
and the quark number scaling of the elliptic flow [10]. 
In fact one of the main criticism to the model was the 
assumption of a quark plasma, which may appear to be 
unlikely due to the preconception of a gluon dominated 
plasma in uRHIC's. 

The impact however of the chemical evolution dis- 
cussed in this Letter appears to be much wider. In fact 
another key physical issue in the field of the QGP physics 
is the determination of the shear viscosity. Recently it 
has been pointed out that a two-component system can- 
not be treated as an effective one-component system with 
an average rj/s from mixture, but one has to solve hydro- 
dynamics for each component [35] . Therefore the relative 
abundance of the two components can significantly af- 
fect the evaluation of the rj/s. Another important issue 
is related with the particle production at high p?- Our 
work shows that a significant gluon to quark conversion 



FIG. 4: Ratio of identified particles from independent frag- 
mentation by AKK parametrizations. The dashed lines are 
the ratio according to the intial quark and gluon distributions, 
the solid lines are after the "chemical" evolution. 



the ratios, (p + p)/2ir°, (p + p)/2K , of identified particle 
can be modified when the gluon conversion mechanism 
is included. This kind of observable asks also for a bet- 
ter knowelcdge of the fragmentation function (FF) in the 
relevant region of parton momentum fraction x and Q 2 . 
Here we have used here the AKK FF [36] just to provide 
a quantitative example of the potential impact of our 
study. This will be however a subject increasing its rel- 
evance with the advent of a new detector for Very High 
Momentum Particle Identification Detector (VHMPID) 
for the ALICE experiment at LHC [37]. 

For high— px physics a further important consequence 
could be the impact on the evaluation of the elastic scat- 
tering energy loss and the gluon radiation mechanism 
responsible for the jet quenching. In fact our study leads 
to a strong decrease of the gg scattering rates and an in- 
crease of the qq and qq that however are known to have 
quite smaller cross sections by a factor 4 — 8 depending 
on the specific channel. Similarly the chemical composi- 
tion of the QGP affects the physics of quarkonia disso- 
ciation in medium. In fact also the scattering rates for 
gJ/fy(T) — s- cc{bb) or g J/^(T) — > gcc(bb) processes can 
be expected to be modified because of the different phase 
space and the reduced scatterer density. The last aspect 
has been indeed very recently investigated in Ref . [38] for 
bottomonium. 

In conclusions our study shows that one could expect 
that the QGP created in uRHIC's, even if it is initially 
a Glasma should very quickly evolve into a plasma dom- 
inated by quark plus antiquarks close to the cross-over 
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temperature T c . The results are quite robust and de- 
volpments of the QPM [19, 20] or inclusion of three-body 
inelastic scatterings may even make the effect larger. The 
result supplies a theorethical justification of the massive- 
quark coalescence hadronization models able to success- 
fully describe several puzzling observations at RHIC and 
LHC; moreover it can have a wide impact on the main 
physical issues of the QGP physics as a quantitative esti- 
mate of the rj/s, the study of the identified high-p^ par- 
ticles and the related jet quenching mechanism as well 
as on the physics of the quarkonia suppression. In gen- 
eral quark dominance along with a small rj/s and a large 



opacity to high-p^ minijets shifts the interpretation of 
the QGP toward even stronger deviations form a per- 
turbative behavior [39] and increased relevance of gluon- 
radiative energy loss [40]. 
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